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This paper reports on a macroscopic fluctuation theory developed over the last ten years 
in collaboration with L. Bertini, A. De Sole, D. Gabrielli and C. Landim. This theory has 
been inspired by and tested on stochastic models of interacting particles (stochastic lattice 
gases). It is the basis for a new approach to the study of stationary non equilibrium states 
applicable to a large class of systems. This overview emphasizes general ideas and for the 
details I refer to the published papers. 

§ 1 . Int r o duct ion 

I will start with a question: why is the theory of irreversible processes so much 
more difficult than the theory of equilibrium phenomena? Here are some reasons. 

In equilibrium we do not have to solve any equation of motion and the Gibbs 
distribution provides the basis for the calculation of macroscopic quantities and their 
fluctuations. In nonequilibrium we cannot bypass the dynamics even in the study 
of stationary states which we may consider as the simplest beyond equilibrium. 
It is easy to write down the equation for the stationary ensemble once we know 
the microscopic dynamics, but the solution is in general very hard to find. Out of 
equilibrium therefore the basic object to define analogs of entropy or thermodynamic 
potentials which we need to construct a natural extension of thermodynamics, is not 
immediately available. 

Since the first attempts to develop a nonequilibrium thermodynamics a guiding 
idea has been that of local equilibrium. This means the following. One assumes 
that locally on the macroscopic scale it is possible to define variables like density, 
temperature, chemical potentials... which vary smoothly on this scale. Microscopi- 
cally this implies that locally the system reaches equilibrium in a time which is short 
compared to the times typical of macroscopic evolutions as described for example by 
hydrodynamics. So what characterizes situations in which this description applies is 
a separation of scales both in space and time. 

A theoretical laboratory where to test this idea is provided by stochastic lattice 
gases which in the last decades have been intensely studied. For several models local 
equilibrium has been proved and hydrodynamic evolution equations have been de- 
rived. The microscopic dynamics of these models is different from the Hamiltonian 
dynamics that one would assume for a classical real gas but the hope is that the 
macroscopic behavior of a system, at least for certain intervals of time, be qualita- 
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tively independent of this assumption. We expect that general properties like the 
type and number of conservations laws should be the most relevant features. 

An important achievement in the study of lattice gases has been the analysis of 
hydrodynamical fluctuations that is the estimate of the probability that the evolution 
of a macroscopic variable, e.g. the density, deviates from a solution of hydrodynam- 
ics. This suggests a dynamical approach to the estimate of the invariant measure in 
analogy with the Freidlin-Ventzell theory for stochastic differential equations in finite 
dimension. Let us give for convenience of the reader a sketch of the Freidlin-Ventzell 
schemed 

Consider a stochastic differential equation 

dxt = b(xt)dt + eadwt, 

where the vector field b is the drift and a the diffusion matrix. One is interested in 
the limit e — > 0. Then the following holds: the probability that the solution stay 
close to a trajectory <pt in a fixed time interval [0,T] is 

P(^~0 t )~exp(-4M^)), (1-2) 

where 

i T (<Pt) = \£ dt{4> - K<t>t))°~ 2 (4> - &(&))• (1-3) 

It is called the large deviation functional. 

From this one obtains the estimate of the stationary distribution in the neigh- 
borhood of an equilibrium point 

P(x)~exp(-iy(x)), (1-4) 

where V{x) = with a trajectory connecting the equilibrium point to x 

and minimizing 1^. If there are several equilibrium points or attractors the theory 
can be easily extended. 

The lesson we learn is that in the limit of small noise an estimate of the stationary 
distribution is reduced to the solution of a variational problem. 

Equation (1.4) reminds of the Einstein theory of equilibrium thermodynamic 
fluctuations^® which states that the probability of a fluctuation from equilibrium 
in a macroscopic region of volume V is proportional to 

exp{VAS/k}, 

where AS is the variation of entropy density calculated along a reversible transfor- 
mation creating the fluctuation and k is the Boltzmann constant. This theory is 
well established and has received a rigorous mathematical formulation in classical 
equilibrium statistical mechanics.^ The analogy with Einstein formula is especially 



Nonequilibrium fluctuations 



3 



suggestive: establishing estimates of this kind in nonequilibrium leads naturally to 
the identification of the exponent with a nonequilibrium thermodynamic function. 
Clearly the role of the small parameter e is taken by the inverse of the volume. This 
is the starting point of the theory reviewed in this paper. For related approaches 
see.™ 

§2. Basic equations and the large deviation functional 

The following equations characterize the macroscopic behaviour of a wide class 
of stochastic lattice gases and agree with the phenomenological equations used so 
far in nonequilibrium thermodynamics near equilibrium, see e.g. the book.^ 

1. The macroscopic state is completely described by the local density p = p(t,x) 
and the associated current j = j(t,x). 

2. The macroscopic evolution is given by the continuity equation 

d t p + V-j = 0, (2-1) 

together with the constitutive equation 

j = J(p) = -D(p)Vp + X (p)E, (2-2) 

where the diffusion coefficient D(p) and the mobility x(p) are d x d positive 
matrices. The transport coefficients D and x satisfy the local Einstein relation 

D(p)=x(p)fo(p), (2-3) 

where /o is the equilibrium free energy of the homogeneous system. 

The equations (|2-1|) - (|2-2|) have to be supplemented by the appropriate boundary 
conditions on dA due to the interaction with the external reservoirs. If \o(x), x € dA, 
is the chemical potential of the external reservoirs, these boundary conditions are 

f!,(p(x)) = X (x), x E dA. (24) 

Some comments are in order. There are no restrictions on the values of the 
chemical potential or on the nonlinearity so that situations far from equilibrium or 
from a stationary state are included. We denote by p = p(x), x € A, the stationary 
solution, assumed to be unique, of equations (|2-ip to (2.4). 

In the constitutive equation (|2-2|) the external field appears linearly which micro- 
scopically means a field of order 1/N acting on each particle, where N is the linear 
dimension of the lattice. A field of order 1 would destroy the diffusive character of 
the evolution giving rise to a hyperbolic evolution equation of first order. 

The transport coefficients D and x depend in general on the density p. In the 
case the system has more than one component, say n, the diffusion coefficient D and 
the mobility become nd x nd matrices. Moreover, in view of Onsager reciprocity, the 
matrix x is symmetric both in the space and in the component indices while D is 
symmetric only in the space indices. In such a case the local Einstein relation (|2-3p 
is D = x R where Rij = d Pi d Pj fo does not depend on the space indices. 
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Summarizing, in the context of stochastic lattice gases, (|2-ip and (|2-2p describe 
the evolution of the density in the diffusive scaling limit, see e.g.®'® The validity of 
the local Einstein relationship (|2-3|) can be deduced from the local detailed balance 
of the underlying microscopic dynamics, see e.g.® 

So far the analysis of the stochastic models can be viewed as a confirmation 
of the phenomenological schemes used by physicists and chemists for a long time. 
The really new information comes from the study of large deviations from the above 
deterministic macroscopic description. 

For a wide class of models the following holds. The stationary measure P s t ad- 
mits a principle of large deviations describing the fluctuations of the thermodynamic 
variable appearing in the hydrodynamic equation. This means the following. The 
probability that in a lattice of N d points the evolution of the so-called empirical 
density p^{Xt) = Y1xgA n Vxi^t) 8 (u— where rj x is the number of particles 
at site x and Xt is the microscopic configuration, deviates from the solution of the 
hydrodynamic equation and is close to some trajectory p(t), is exponentially small 
and of the form 

P st ( PN (X t ) ~ p(t),te [hM) ~ e-^FWW+'toMW, (2-5) 

where I(p) is a functional which vanishes if p(t) is a solution of (|2-ip - (j2-2p and 
J-~(/5(ii)) is the cost to produce the initial density profile p(ti). We normalize it so 
that J-{p) = 0. The functional I{p) represents the extra cost necessary to follow 
the trajectory p{t). Finally p^{Xt) ~ p{t) means closeness in some metric and ~ 
denotes logarithmic equivalence as N — > oo. 

A simple argument based on time reversaP® now shows that 

J-( /0 ) = inf/ hoo . 0] ( / 5), (2-6) 

p 

where the inf has to be taken over all trajectories connecting the stationary state p 
to the profile p. The same argument identifies the minimizing trajectory, see later. 
The quantity J-(p) is an infinite dimensional analog of the V{x) in Freidlin-Ventzell 
theory and has to be identified with a non equilibrium thermodynamic function. The 
fluctuations we are considering take place at constant temperature and volume so 
that it is reasonable to identify J-(p) with the variation of the nonequilibrium free 
energy to produce the profile p. If the stationary state is an equilibrium state the 
Einstein theory is recovered as it will be discussed later. 

The explicit formula of the functional Jr tlit2 i(p)] is 

I[ Tl ,T 2 ] (P) = \ J^dt ( [d t p + V • J(p)] Kip)- 1 [d t p + V • J(p)] ), (2-7) 

where the positive operator K(p), the analog of a 2 in Freidlin-Ventzell theory, is de- 
fined on functions u vanishing at the boundary dA by K(p)u = —V • (x(p) Vit) . The 
expression (2.7) is similar to (1.3), i.e. the fluctuations of the thermodynamic vari- 
ables of a stochastic lattice gas are formally those of a stochastic partial differential 
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equation where the deterministic part is given by hydrodynamics. These equations 
are known as fluctuating hydrodynamics and are often used in phenomenological 
calculations. One has to be careful however because it is often difficult to attribute a 
precise mathematical meaning to stochastic nonlinear partial differential equations. 
Divergences akin to those in quantum field theory appear in dimension greater than 
1. 

For a heuristic derivation of (2.7) seeP^ There is factor 1/2 of difference due to 
a slightly different definition of the transport coefficients D and \- 

§3. The Hamilton-Jacobi equation and the interpretation of T 

As shown iiPSJffi) the functional J-, as defined in (2.6), is the maximal solution 
of the infinite dimensional Hamilton-Jacobi equation 

(Vf -*)vf)-(fv., W )^0, (3,, 

where, for p that satisfies ()2-4f) . 5J~/5p vanishes at the boundary of A. At the 
macroscopic level this condition reflects the fact that we consider variations of the 
density that do not change the boundary values. The arbitrary additive constant on 
the maximal solution of (|3-1|) is determined by the condition J-(p) = 0. By maximal 
solution we mean that any solution F(p) to (|3-ip (satisfying F{p) = 0) is a lower 
bound for T . In nonequilibrium we expect generically T to depend nonlocally on 
the density. 

It is easy to see that the Hamilton-Jacobi equation implies that VJ can be 
decomposed in such a way that the hydrodynamic equations can be written 



d t p = -VJ(p) = V • (x(p)V— ) + A(p), (3-2) 

where A(p) satisfies the orthogonality condition 

(¥-,A(p)) = 0. (3-3) 

This decomposition confirms the interpretation of J 7 as a nonequilibrium free 
energy. In fact the first term is the dissipative part of the equation which determines 
the relaxation to the stationary state due to the thermodynamic force V^p while 
A(p) does not contribute to the spontaneous relaxation. A(p) is minus the divergence 
of the non-dissipative part of the current and vanishes in the stationary state. It 
can be shown thatf^ the hydrodynamics associated to the time reversed microscopic 
process can be written 



d t p = -vr(p) = v • (x(p)v— ) - A( P ), 



(3-4) 
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where J*{p) is the current associated to the time reversed process. By summing 
equations (3.2) and (3.4) we obtain the nonequilibrium fluctuation-dissipation rela- 
tion 

V(J*(p) + J(p)) = -2V X (p)V 6 -f. (3-5) 

dp 

Actually a similar relation holds for the currents^} 

r(p) + J(p) = -2 X (p)V 6 -f. (3-6) 

dp 

We can now identify the optimal (minimizing) trajectory in (2.6). Let J- be the 
maximal solution of the Hamilton- Jacobi equation and J* as in (|3-6p . Fix a time 
interval [Ti,T2] and a path p(t), t S [Ti,T2]. We claim that 

I [TljT2] (p)=T{p(T 2 ))-F(p(T 1 )) 

+\JJdt ([dtp-V ■ J*(p)] K{p)- X [dtp-V ■ J*(p))), (3-7) 

as can be shown by a direct computation using (2.7), the Hamilton- Jacobi equation 
(|3-ip and (3.6). From the identity (|3-7p we immediately deduce that the optimal 
path for the variational problem (2.6) is the time reversal of the solution to (3.4) 
with initial condition p. 

We next show that according to definition (2.6), the free energy J-(p) is equal 
to the work done by the thermodynamic force on the system along the optimal 
trajectory p*. Indeed, by the above identification of p* 

dt (jp-' dtp * 

dtff iV .jV)}= f dt((-r( P *))-v 5 -f). 
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The decompositions (3.2)-(3.4) remind of the electrical conduction in presence 
of a magnetic field P*^ Consider the motion of electrons in a conductor: a simple 
model is given by the effective equation,^ 

p = -e(E + — p AH) - -p, (3-8) 
V mc / r 

where p is the momentum, e the electron charge, E the electric field, H the magnetic 
field, m the mass, c the velocity of the light, and r the relaxation time. The dissi- 
pative term p/r is orthogonal to the Lorenz force p A H. We define time reversal as 
the transformation p i— > — p, H i-> — H. The time reversed evolution is given by 

p = e(E+ — p AH) - -p. (3-9) 
V mc J t 

Let us consider in particular the Hall effect where we have conduction along a rectan- 
gular plate immersed in a perpendicular magnetic field H with a potential difference 
across the long side. The magnetic field determines a potential difference across the 
short side of the plate. In our setting on the contrary it is the difference in chemical 
potentials at the boundaries that introduces in the equations a non-dissipative term. 



Nonequilibrium fluctuations 



7 



§4. Equilibrium states and macroscopic reversibility 

In this section we consider the notion of equilibrium from the standpoint of 
nonequilibrium^^ We define the system to be in equilibrium if and only if the 
current in the stationary profile p vanishes, i.e. J{p) = 0. Nonetheless, in presence 
of external (electric) fields and boundary reservoirs an equilibrium state can be 
highly inhomogeneous. An example of such a situation is provided by sedimentation 
equilibrium in gravitational and centrifugal fields. In spite of this, the free energy is 
a local function of the state variables and coincides locally with the equilibrium free 
energy in absence of external fields and boundary driving. In particular there are no 
macroscopic correlations. 

For an equilibrium state, characterized by a constant or non constant stationary 
profile p(x) such that J(p) = 0, the free energy functional J- is obtained as follows. 
Let 

f(p,x)= f dr f dr'f(j(r') = f (p)-f (p(x))-f( ) (p(x))[p-p(x)}, 

Jp{x) Jp(x) 

where, we recall, /o(p) is the equilibrium free energy density. It is easy to shovP* 
that the maximal solution of the Hamilton-Jacobi equation is 

F{p)= [ dxf(p(x),x). (4-1) 

J A 

We emphasize that the above argument depends crucially on the structure (|2-2p of 
the current and on validity of the local Einstein relation (|2-3p . 

The condition J(p) = can be showrPS to be equivalent to J*(p) = J(p) for 
any profile p. We call the latter condition macroscopic reversibility. The notion 
of macroscopic reversibility does not imply that an underlying microscopic model 
satisfies the detailed balance condition. Indeed, as it has been shown by explicit ex- 
amples''^ 1 there are non reversible microscopic models which are macroscopically 
reversible. If the system is macroscopically reversible the optimal trajectory to create 
a fluctuation is the time reversal of the relaxation trajectory solution of (|2-l|) - (|2-2p . 
IrP we called this symmetry the Onsager-Machlup time reversal symmetry.®'^ 

We have defined the macroscopic reversibility as the identity between the cur- 
rents J{p) and J*{p)- We emphasize that this is not equivalent to the identity 
between V • J(p) and V • J*(p). Indeed, it is possible to give examples 151 of a non 
reversible system, i.e. with J(p) ^ 0, such that the optimal trajectory for the varia- 
tional problem (2.6) is the time reversal of the solution to the hydrodynamic equation 

Our analysis of equilibrium so far does not include magnetic fields but can be 
extended to cover this situation. It can be done by distinguishing dissipative and 
non dissipative currents: a natural definition of equilibrium is then the vanishing of 
the dissipative currents. 
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§5. Long range correlations 

Space correlations extending over macroscopic distances appear to be a generic 
feature of stationary nonequilibrium states, a fact known since a longtime.®® The 
ultimate reason for these correlations seems to be connected with the violation of time 
reversal as they can appear also in equilibrium states if the microscopic dynamics 
strongly violates time reversal invariance.^'^J An interesting question is to what 
extent long range correlations make nonequilibrium akin to critical phenomena. 

In our approach, since we are not limited to the vicinity of the stationary state 
we can obtain the (nonequilibrium) density correlations functions of arbitrary order 
in terms of the functional derivatives of T . In general the functional T cannot be 
written in a closed form, but - by a suitable perturbation theory on the Hamilton- 
Jacobi equation (|3- 1[) - we can derive such correlations functions. In this section 
we discuss the two-point correlation and establish a criterion to decide whether the 
density fluctuations are positively or negatively correlated. Recursive equations for 
the correlation functions of any order are established inP^ We emphasize that we 
are concerned only with macroscopic correlations which are a generic feature of 
nonequilibrium models. Microscopic correlations which decay exponentially or as a 
summable power law disappear at the macroscopic level. 

We introduce the pressure functional as the Legendre transform of free energy 

T 

G(h) = sup {(hp) -Hp)}- 
p 

By Legendre duality we have the change of variable formula h = p = so 
that the Hamilton-Jacobi equation (|3-ip can be rewritten in terms of Q as 



(V,.,(|)V A )-(V,.,(|)V|-,(|) £ ).0, <M) 

where h vanishes at the boundary of A. As for equilibrium systems, Q is the gener- 
ating functional of the correlation functions. In particular, by defining 



c(,,»)- s " g{h) 



Sh(x)5h{y) 
we have, since J- has a minimum at p, 



h=0 



g(h) = (h,p) + ±(h,Ch) + (h 2 ), 



Hp) = \{(p~ Pi c-\ p - p)) + o(( P - pf). 



or equivalently 



By expanding the Hamilton-Jacobi equation (|5-ip to the second order in h, and 
using that 5Q /5h(x) = p(x) + Ch(x) + o(h 2 ), we get the following equation for C 



V/i- 



X (p)Vh - V{D{p)Ch) + x'(p)(Ch)E ) = 0. (5-2 
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We now make the change of variable 

C(x, y) = C cq (x)5(x -y) + B(x, y), 

where C eq (x) is the equilibrium covariance. By using (|2-3p we deduce that 

C eq {x) = D- l {-p{x))x(p{x)). 

Equation (|5-2|) for the correlation function then gives the following equation for B 

C ] B(x,y) = a(x)5(x-y), (5-3) 

where & is the formal adjoint of the elliptic operator C = L x + L y given by, using 
the usual convention that repeated indices are summed, 

L x = D tJ (p(x))d Xi d X] + ^{p{x))E 3 {x)d Xi , (5-4) 

and 

= d Xi MM*)) D^{p(x))J k (x)] , 

where we recall J = J(p) = —D(p(x))Vp(x) + x(p( x ))B(x) is the macroscopic 
current in the stationary profile. 

In equilibrium systems J = so that we have a = 0, hence B = 0, namely there 
are no long range correlations and C(x,y) = C eq (x)5(x — y). Moreover, since C is 
an elliptic operator (i.e. it has a negative kernel), the sign of B is determined by 
the sign of a: if a(x) > 0, Vx, then B(x,y) < 0, Vx,y, while if a(x) < 0, Vx, then 
B(x,y) > 0, Vx,y. For example, consider the following special case. The system 
is one-dimensional, d = 1, the diffusion coefficient is constant, i.e. D{p) = Dq, the 
mobility \{p) ls a quadratic function of p, and there is no external field, E = 0. 
Then 

B(x, y) = -^- X "{yp) 2 ^\x, y), (5-5) 

where A~ 1 (x,y) is the Green function of the Dirichlet Laplacian. Two well studied 
models, the symmetric exclusion process, where x(p) = p(l — P)i an d the KMP 
process,^ 1 where x(p) = P 2 , meet the above conditions. Then (|5-5p shows that their 
correlations have opposite signs. 

§6. Thermodynamics of currents 

In nonequilibrium a very important observable is the current flux. This quantity 
gives information that cannot be recovered from the density because from a density 
trajectory we can determine the current trajectory only up to a divergence free vector 
field. To discuss the current fluctuations, we introduce a vector-valued observable 
J7at({^(t, < a < t}) of the microscopic trajectory Xt which measures the local net 
flow of particles and satisfies formally the microscopic continuity equation 



dtPN + Vat • J N = 0, 
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where Vat is the gradient on the lattice. For the details seeP^'® As in the case 
of the density, for stochastic lattice gases, we are able to derive a dynamical large 
deviations principle for the current. If Px stands for the probability P s t conditioned 
on the initial microscopic state Xq, given a vector field j : [0, T] x A — >• M. d , we have 

P Xo [Jn{X) « j(t,u)) ~ exp{ - N d l [0>T] (j)}, (6-1) 

where the rate functional is 

%T](j) = \£dt([j- Jip^xipr'U-Ap)]). (6-2) 

We recall that 

J(p) = -D(p)Vp + X (p)E . 

In (6.2) p = p(t, u) is the solution of the continuity equation dtp + V • j = with 
the initial condition p(0) = po associated to Xo. The rate functional vanishes if 
3=J(P)- 

Among the many problems we can discuss within this theory, the fluctuations 
of the time average of the current Jn over a large time interval have been analysed. 
This question was addressed by Bodineau and Derrida irPSJ in one space dimension 
by postulating an "additivity principle" which relates the fluctuation of the time 
averaged current in the whole system to the fluctuations in subsystems. However 
their approach does not always apply. In fact the probability of observing a given 
divergence free time averaged fluctuation J can be described by a functional 
which we characterize, in any dimension, in terms of a variational problem for the 
functional 2jo,r] 

<P(J) = lim infiz [0 , T] (i), (6-3) 

where the infimum is carried over all paths j = j(t,u) having time average J. The 
static additivity principle postulated in® gives the correct answer only under addi- 
tional hypotheses which are not always satisfied. Let us denote by U the functional 
obtained by restricting the infimum in (|6-3|) to divergence free current paths j, i.e. 

U{J) = ind([J- J{p)], X {p)- l [J - J(p)}), (6-4) 
P 4 

where the infimum is carried out over all the density profiles p = p(u) satisfying the 
appropriate boundary conditions. From (|6-3p and (|6-4p it follows that <P < U. In 
one space dimension the functional U is the one introduced inP^ While <P is always 
convex the functional U may be non convex. In such a case U(J) underestimates 
the probability of the fluctuation J. IrJI3'l2I} W e interpreted the lack of convexity of 
U, and more generally the strict inequality ^ < U, as a dynamical phase transition. 

There are cases in which <E> = U. Sufficient conditions on the transport coef- 
ficients D, x f° r the coincidence of <P and U can be givenP^J Consider the case 
when the matrices D(p) and x{p) are multiples of the identity, i.e., there are strictly 
positive scalar functions still denoted by D(p), x(p)> so that D(p)ij = D(p)5ij, 
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x(p)i,j = x{p)&i,j-> i, j = 1, ■ ■ ■ ,d. Let us first consider the case with no external field, 
i.e. E = 0; if 

D(p) X "(p) < D'(p) X '(p), for any p, (6-5) 

where ' denotes the derivative, then ^ = U. In this case U is necessarily convex. 
Moreover if 

D {p)x"{p) = D'(p)x'(p), for any p, (6-6) 

then we have <P = U for any external field E. 

To exemplify situations in which <P < U consider the fluctuations of the time 
averaged current for periodic boundary conditions. Two models have been dis- 
cussed so far. The Kipnis-Marchioro-Presutti (KMP) model^S which is defined 
by a harmonic chain with random exchange of energy between neighboring oscilla- 
tors, and the exclusion process. In the case of the KMP model we have U(J) = 
(l/4)J 2 /x(m) = (l/4)J 2 /m 2 , where m is the (conserved) total energy. For J large 
enough, ^(J) < U(J). This inequality is obtained by constructing a suitable trav- 
elling wave current path whose cost is less than f/(J)P2} For J not too large the 
additivity principle holds as it has been verified numerically in.® A similar result 
has been obtained by Bodineau and DerridgP^ 1 for the periodic simple exclusion pro- 
cess with external field. For the KMP process this phenomenon is rather striking as 
it occurs even in equilibrium, i.e. without external field. 

The behavior of I and & under time reversal shows that satisfies a fluctuation 
relationship akin to the Gallavotti-Cohen theorem for the entropy production.® - ^} 
The anti-symmetric part of tf> is equal to the power produced by the external field 
and the reservoirs independently of the details of the model 

#(J) -<£>(- J) =#(J) -$ a {J) = -2{J,E) + / dUXoJ-h, (6-7) 

JdA 

the right hand side of this equation is the power produced by the external field and 
the boundary reservoirs (recall E is the external field and Ao the chemical potential 
of the boundary reservoirs). From this relationship one derives a macroscopic version 
of the fluctuation theorem for the entropy production. 

For recent interesting results obtained from the macroscopic fluctuation theory 
in the study of current fluctuations seeP^'GS 

§7. Conclusions and comparison with other approaches 

The theory developed so far, as emphasized inpJ can be viewed as a selfcon- 
tained macroscopic description of diffusive systems out of equilibrium. It allows a 
clear identification of the dissipative part and of thermodynamic forces in relaxation 
phenomena. Long range space correlations among thermodynamic variables are a 
generic consequence of the theory and equations for the correlation functions of any 
order have been established.^ Also different regimes in current fluctuations are 
predicted that have been characterized as dynamical phase transitions. These tran- 
sitions are shown to exist in the simplest models considered and it is a challenge to 
discover them in real systems. 
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The thermodynamic functionals J- and ^ can be calculated from dynamical 
measurable quantities like the transport coefficients D and x- This is a deep differ- 
ence with respect to equilibrium thermodynamics where static properties like specific 
heats are involved. Prom the standpoint of nonequilibrium it is natural to consider 
also the equilibrium free energy as determined by the transport coefficients through 
the Einstein relation. 

Another substantial difference is that in the definition of nonequilibrium ther- 
modynamic functionals optimal trajectories are involved which are different from the 
infinitely slow reversible transformations of classical thermodynamics: the optimal 
trajectories go through nonequilibrium nonstationary states. In the theory developed 
so far the boundary conditions are kept fixed: the study under boundary conditions 
(chemical potentials, volume...) which slowly vary on the macroscopic time scale is 
a next natural step. 

A different very general approach to the theory of nonequilibrium stationary 
states was initiated by Oono and PaniconP^ anc j pursued in the work of Hatano, 
Hayashi, Sasa, Tasaki.® - ® while ijJSIiEB microscopic models underlie the analysis, 
irP« a guiding idea is to keep as much as possible the phenomenological character of 
classical thermodynamics without reference to an underlying microscopic dynamics. 
In this work the authors discuss the operational definition of nonequilibrium ther- 
modynamic observables in concrete situations and generalize basic operations like 
decomposition, combination and scaling of equilibrium thermodynamics to nonequi- 
librium states. The possibility of experimental tests is then discussed. 

In more recent papers by Komatsu, NakagawsP^ 1 and by Komatsu, Nakagawa, 
Sasa, TasakiJSSJfiU the problem of constructing microscopic ensembles describing 
stationary states of both stochastic and Hamiltonian systems, is considered. Expres- 
sions for the nonequilibrium distribution function are proposed either exact or valid 
up to a certain order in the parameters keeping the system out of equilibrium. The 
key quantity appearing in these expressions is the entropy production. In® heat 
conduction is considered in particular and an extension of the Clausius and Gibbs 
relations is derived. 

A direct comparison between these works and ours is not immediately available 
as we are asking different questions but an effort should be made to construct a 
bridge between them. 

In a recent paper Ottinger® compared his own approach® to nonequilibrium, 
called GENERIC (general equation for the nonequilibrium reversible-irreversible cou- 
pling), with our macroscopic fluctuation theory discussing aspects where a correspon- 
dence could be established. His starting point is a separation in the macroscopic 
evolution equations of dissipative and conservative terms which reminds of our de- 
composition (3.2). An important difference is related to the fact that our free energy 
out of equilibrium is generically nonlocal in space and this is connected with the ex- 
istence of long range correlations. There is more to understand and the comparison 
should be developed further. 
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